Coherent molecular solitons in Bose-Einstein condensates 
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We analyze the coherent formation of molecular BEC from 
an atomic BEC, using a parametric field theory approach. We 
point out the transition between a quantum soliton regime, 
where atoms couple in a local way - to a classical soliton 
domain, where a stable coupled-condensate soliton can form 
in three dimensions. This gives the possibility of an intense, 
stable atom-laser output. 
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Parametric sohtons or simultaneous solitary waves 
( "simultons" ) , involving the optical x*"^"* nonlinearity, 
have been the topic of much recent theoretical and exper- 
imental interest in nonlinear optics. We propose a novel 
mechanism by which a similar phenomenon may occur 
in nonlinear atomic optics, in which coherent molecule 
formation in a Bose-Einstein condensate takes the place 
of second harmonic generation. 

This requires a coupling that converts two atoms into 
one molecule, thus generating coupled atomic and molec- 
ular Bose-Einstein condensates - and so taking advantage 
of molecular states that are known to exist in alkali-metal 
vapors. Our model includes a coherent molecular forma- 
tion process (i.e., without dissipation) in an atomic BEC 
vapor (or atom laser [||), either through a Feshbach 
resonance Q or Raman photo-association We note 
that Feshbach resonances have already been observed . 
The coherently coupled atom-molecular condensate could 
provide a route to the observation of a localized three- 
dimensional BEC soliton, even in the absence of a trap 
potential. A possible application is in the free propaga- 
tion of a non-diverging atom laser pulse, thus greatly in- 
creasing the intensity in an atom laser beam. Even more 
than this, would be the importance of observing the strik- 
ing physical properties of this novel quantum field theory, 
and the corresponding Bose-enhanced chemical kinetics. 

The original solution for the parametric soliton was 
in a one-dimensional environment These classical 
solutions have been classified topologically 0, and are 
generic to the mean-field theories of parametric non- 
linearities that convert one particle into two (and vice- 
versa). The equations are non-integrable, and are differ- 
ent to the usual integrable classes of soliton equations. A 
considerable advantage of these types of nonlinear equa- 



tions is that they are capable of providing solutions in 
one, two, or three space dimensions, which does not oc- 
cur in the usual Gross-Pitaevskii equations. Both clas- 
sical [^-^ and quantum |^ solutions have been recently 
identified (including observation of classical solitons in 
experiment 0), although these different types of soli- 
ton have strikingly different qualitative behavior. 

The purpose of this Letter is to point out the physical 
origin of these differences between the quantum and clas- 
sical versions of the parametric field theory; and to iden- 
tify experimental requirements for observing these novel 
effects in Bose condensates. We consider the following 
basic Hamiltonian, to give a simple model of molecule 
formation: 



H — Ho + Hi + Hint , 
where the free and interacting Hamiltonians are: 
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Here we define complex fields $ = / d^ka(k) exp[i(k- x)] 
and * = / (i^k6(k) exp[i(k • x)] . The field $ represents 
an atomic species of mass m in a potential V$ (x) , in one 
internal state, while ^ represents a dimer species of mass 
M = 2m, in a single vibrational and rotational state, 
with a potential V3,(x). 

The coupling constant x represents a formation rate 
for the dimer, in the S-wave scattering limit, while k rep- 
resents the effective self- interaction of the atomic field. In 
the absence of any trap, the potentials are uniform, and 
hp = h{V\i/ — 2V$) is the formation energy of the dimer 
species. We note that these interactions are idealized, in 
the sense that both x s^nd k represent processes that are 
microscopically nonlocal. To represent such nonlocal be- 
havior, we must introduce a momentum cutoff km in the 
relative momenta of interacting fields, which physically 
must be around the inverse S'-wave scattering length - if 
we wish to use the non-renormalized effective potential 
to describe S'-wave scattering. This is known to be essen- 
tial to the correct interpretation of these types of effec- 
tive field theories. It should be recognized that molecular 
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self-interactions - as well as atom-molecular scattering - 
will occur as well. These are neglected here, since the 
relevant cross-sections are not well known. 

In the corresponding nonlinear optical case, the <i> and 
^' fields would correspond to a first and second harmonic, 
coupled by a x^^^ nonlinearity of the dielectric, while n 
would correspond to a x'-'^-' nonlinearity. The interplay 
between quadratic and cubic nonlinearities in the case 
of nonlinear optical solitons has been analyzed, at the 
classical level and for one space dimension, in The 
effective masses, which should be different in the longi- 
tudinal and transverse directions, describe the effects of 
dispersion and diffraction, respectively, for both the fields 
(see, e.g., § for more details). Here the equations refer 
to a moving frame situation, with coordinates moving at 
the group velocity. 

By comparison, in the directly comparable atomic case, 
we are considering atoms in free-space. No potential 
needs to be included, since this is not essential to soliton 
formation. The molecular formation process would be 
tuned in any practical experiment, by magnetic fields or 
external Raman coupling, in order to reduce the energy 
mismatch hp. An important consideration is the possible 
effects of losses due to inelastic atom-molecule collisions. 
We assume that an appropriate choice of molecular levels 
is made, so that these losses can be ignored over the rele- 
vant time-scales for solitons to form. Thus, the neglect of 
molecular vibrational transitions is crucial to the present 
theory, which only includes one molecular level. An ideal 
situation would involve a direct coupling via a tuned Ra- 
man transition to the molecular ground state. A more 
sophisticated theory would include detailed atomic po- 
sitions and multiple energy levels within each molecule. 
Our theory neglects these additional complications. 

The Heisenberg equations of motion that correspond 
to the basic Hamiltonian are: 

^-V2$ + x*$^ + + K(,(x)l> , 
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As a first step, we can take mean values, so that 
(j> — {<t) and = ($), and assume operator product 
factorization. This gives rise to mean field equations, 
valid for a momentum cutoff less than the S-wave scat- 
tering length. For the case of Bose condensates in exist- 
ing evaporative cooling experiments, near the atomic col- 
lective ground state, the mean field equations represent 
modified Gross-Pitaevskii equations - which are known 
to successfully describe BEG excitations. 

Another way to understand the behavior of this quan- 
tum many-body system is to look for energy eigenstates 
of the original Hamiltonian, in the limit of a large mo- 
mentum cutoff. These must simultaneously be the eigen- 



states of N = J d^x 
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conserving the gen- 



count each molecule as two atoms). Solving this, a re- 
markable fact emerges. We can show rigorously that in 
the limit of free space propagation, an A^-boson ground 
state exists - by finding exact upper and lower bounds 
on the Hamiltonian energy. Since these coincide in three 
dimensions, we have the result that the (idealized) quan- 
tum ground state energy is exactly: 



E. 



N 



N ( hx^ 
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where we assume N is even. The proof of the lower bound 
also assumes k > and x^ > 2pK, and the result is 
obtained using the known solution of the two-particle 
(N = 2) bound-state problem 

This corresponds to N/2 independent quantum soli- 
tons or "dressed" molecules, each of which exist in a lin- 
ear superposition with a pair of atoms (like a Cooper 
pair), so that: 



N\ 
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d'^k g(k)aT(k)aT(-k) 
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In this limit of a large cutoff in the quantum field theory, 
the ground-state energy has no lower bound as k —^ 0. 
This is in remarkable contrast to the known mean-field 
behavior of the corresponding classical energy, which is 
rigorously bounded below (see, e.g., [^). Of more inter- 
est is the limiting behavior of the ground-state quantum 
energy when there is a cutoff fcm present. We have ob- 
tained a variational estimate of this quantity, and for this 
case we obtain (p, k — > 0): 



(6) 



Here we have taken the case of a relatively large cutoff, 
so that the result assumes that km ^ {xm / {2TTh))'^ , and 
uses a variational ansatz of the form given previously. 
The ansatz gives us the true ground state energy in the 
limit km oo, (for any finite k), since upper and lower 
energy bounds coincide. However, it is not necessarily 
the lowest possible energy at finite km- In order to show 
this, we consider a coherent or mean-field ansatz, with 
broken symmetry, of form: 



exp 



(x)$t(x)+^(x)*t(x) l|0) 
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eralized particle number N (total number of atoms if we 



For this case, the classical decorrelation originates in 
coherent-state factorization properties of the Hamilto- 
nian. This state is, however, not an eigenstate of H 
(since it is not an eigenstate of N). It is an approximate 
(semi-classical) eigenstate at large N, and corresponds 
to two coupled Bose-Einstein condensates under broken 
symmetry conditions. 
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We will now show that, provided tp{yi), (/'(x) are cho- 
sen to minimize the classical Hamiltonian, they can give a 
lower energy than previously - although still finite. This 
calculation makes use of the known result that the classi- 
cal parametric Hamiltonian is always bounded below , 
and the bound is given by the soliton energy for exact 
phase matching p = 0. This soliton energy is estimated 
by means of a variational ansatz applied to the Hamilto- 
nian. We choose 

0(x) - .giiV2[2/(7rsi)]3/4exp(-|xp7VVsi) , 

V'(x) - -g2N'[2/{7:s2)?^^ exp(-|xpiVVs2) . (8) 

The negative sign for '(/'(x) ensures that the coupling 
energy is negative, and the normalization implies that 
51+52 — 1- We note that although a uniform varia- 
tional ansatz is possible, it is known that a uniform field 
of this type is always unstable for a purely parametric 
coupling - and hence cannot give the lowest energy. 
Substituting into the Hamiltonian, gives us the result: 
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where we have used M = 2m, with the simplified no- 
tation of X = 2^/2(2/7r)3/'*mx/(3?i) and, similarly: k = 

2-5/2(2A)3/2k.^ 

To minimize i?^, under the constraint of a fixed N ^ is 
a non-trivial algebraic procedure. However, the physics 
is considerably simplified in the region where the term in 
iV'^ is dominant - which we note should not involve too 
large a contribution from the repulsive term that scales 
with A^^, and tends to destabilize soliton formation. In 
this region (i.e., assuming k ~ p ~ 0), we obtain a cou- 
pled molecular Bose condensate minimum energy of: 



This question is therefore a subtle combination of mo- 
mentum cutoff and particle density effects. To give some 
numerical results we consider m ~ lO^^sj^g^ a,nd use a 
X- value estimate of about x ~ 10~^ m'^/^/ggc (given in 
1^, by Tommasini ei. al., for a Feshbach resonance Q), 
leading to mx/Ti ~ 10'^ m^^/2_ With a choice of the cut- 
off at km ^ lnm~^, this gives a critical atom number 
of about Ncr ~ 10^, which is well within the range of 
current BEC experiments. At low particle density, the 
formation of individual dressed molecules is favored, as 
atoms couple to molecules in a particle-like way. The pro- 
cess is analogous to Rabi oscillations of atoms between 
two different electron sub-levels, except that is occurs 
between pairs of atoms and the corresponding molecular 
levels. These dressed states have interesting properties, 
reminiscent of Cooper pairs, but cannot be described by 
the classical parametric soliton equations. At large cou- 
plings X, and at large density (but not too large so that 
5-wave scattering is dominant) the coherent coupling of 
two entire condensates is dominant - just as in nonlin- 
ear optics. In this domain, provided other recombination 
processes are negligible, there are strong, coherent and 
nonlinear wave-like interactions between the atomic and 
the molecular Bose condensates. For these parameters, it 
even appears possible to form a stable, three-dimensional, 
Bose-Einstein soliton. 
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where C is a constant given by C ~ 1 • 2 x 10"^. The rel- 
evant length scale is nearly identical for the two coupled 
condensates, and is given by: 
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This enables us to more clearly understand the appar- 
ent paradox that a full quantum theory gives a qualita- 
tively different lower energy bound to the corresponding 
classical mean field theory. To obtain a stable coupled 
atom-molecular condensate, we require E^ < Eq , which 
occurs at a critical boson number: 



N > Nr^ 



(12) 
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